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We develop a general analysis of the diffusive dynamics of polydisperse polymers in the presence 
of chemical potential gradients, within the context of the tube model (with all species entangled). 
We obtain a set of coupled dynamical equations for the time evolution of the polymeric densities 
in a form proposed phenomenologically in recent work by N. Clarke, but with explicitly derived 
coefficients. For the case of chemical polydispersity (a set of chains that are identical except for 
having a continuous spectrum of enthalpic interaction strengths) the coupled equations can be fully 
solved in certain cases. For these we study the linearised mode spectrum following a quench through 
the spinodal, with and without a passive (polymeric) solvent. We also study the more conventional 
case of length polydisperse chains in a poor solvent. Here the mode structure is more complicated 
and exact analysis difficult, but enough progress can still be made to gain some qualitative insight, 
ryj ' We briefly discuss the modifications required to allow for the presence of unentangled, low molecular 

weight species in the system. 

I. INTRODUCTION 

There has been significant recent progress in the study of polydisperse polymers under equilibrium conditions. 
£j ' This includes work on static phase coexistence (reviewed in ref 1) and on interfacial thermodynamics in situations 
of spatial inhomogeneity ||. Much of this work has been based on exploiting the fact that, within a standard 

■ Flory Huggins picture of the polymeric systems involved, the excess free energy depends only on a small number of 
moments of the distribution of polymer lengths or, if the chains are polydisperse not in length but in some other 
feature, moments of the distribution of that feature. An example of the latter is chemical polydispersity, where chains 

I/"") | have, in effect, a continuous distribution of \ parameters Q but are otherwise identical; see also ref 3. (This can be 
. realized experimentally by having a random copolymer in which the fraction of each type of monomer varies from 
chain to chain.) This moment structure carries over into the nonlocal excess free energy for inhomogeneous states, at 
least in limiting cases Q ; the moments become local moment density fields which depend nonlocally on the number 
densities of the individual polymeric species. 

Much less is understood about the dynamics of polydisperse systems. Consider, for example, the early stages of 

■ spinodal decomposition following a quench into the unstable regime. For length polydispersity this has been studied 
in recent papers by Clarke Q and by Warren [^), but even the qualitative picture remains far from clear. One 

I . complicating factor behind this is that, for the normal case of length polydispersity, the simplifications brought about 
by the moment structure in the excess free energy do not lead to any obvious simplification of the dynamical equations: 
the natural coordinates for dynamics are not the moment densities, but the species monomer densities. (See ref 6 for 
a related discussion, in the context of colloids, addressed there via a perturbation theory in the narrowness of the size 
distribution.) Previous treatments of length polydispersity mapped the polydisperse situation into that of a mixture 
with a finite set of components by matching the lowest moments of the distribution [Q . For chemical polydispersity, 
however, there is significantly more carry-over of the simplifications in the excess free energy to dynamics. The basic 
reason for this is that such chains are all dynamically degenerate (they have the same mobilities and structures) even 
though they are subject to different thermodynamic forces. Hence they can, in effect, be linearly recombined into 
moment densities, which diagonalize the thermodynamics. This is not possible for length polydispersity where the 
mobilities (and also the structure factors) of the polymeric species vary from one chain to another; in general there 
is then no practical way to simultaneously diagonalize the thermodynamic and the mobility-related factors in the 
diffusion matrix. 

In this paper we develop a general formalism which includes, in principle, both chemical and length polydisper- 
sity (among others). In part this formalism resembles that of Clarke J4| but it is more general and supplements 
his phenomenological equations with explicit calculation of the coefficients from a microscopic model based on the 
tube dynamics. The latter approach was developed previously in a similar context (but without polydispersity) by 
Brochard ||. These developments are presented in sections jf] and III. In section IV we present the free energy 
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functional employed thereafter and in section |y| we address in some detail the case of chemically polydisperse chains. 
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These are considered in the melt and also in the presence of a passive solvent (one that does not have any direct 
enthalpic preference between the chemically polydisperse species). For the first case the dynamical equations can be 
solved in detail; in the second case, we find a relatively simple mode structure requiring a two-by-two diagonalization 
whose details involve the shape of the polydisperse distribution in the homogeneous (parent) phase. This simplifies 
further for the case of a symmetric parent. 

These results are then contrasted with those in the literature for the case of length polydispersity, which is studied 
further in section [VT| . Although, for the reasons stated above, the same level of progress cannot be made in this 
case, the formal structure of the problem (which involves, for example, a coupling to negative moments of the length 
distribution) is rather interesting. We are able to make some indicative progress by usin g a truncation scheme along 



the lines suggested by Warren H, and indicate avenues for further analysis. Section VII contains our conclusions 



II. SYSTEM AND BASIC DYNAMICS 

We consider a polydisperse mixture of polymers. Each species is characterized by the chain number density pi(r). 
(We use a discrete notation, but the continuum limit can be taken finally, as required.) The corresponding monomer 
density is written as 4>i(r) = NiPi(r), where iVj is the polymerization index of species i; note that this neglects terms 



0(V ) terms, see sec.1V. Depending on the kind of polydispersity, Ni may be the same for different species. When 
dealing with polydisperse mixtures, the chain number density is the natural density to work with (as opposed to 
monomer density), because it characterizes the overall motion of the chain. 
The continuity equation for polymer species i is 

| = -V-(p iVl ) (1) 

where v, is the mean velocity of species i, and pjVj = J' is the total flux of that species. We can rewrite this equation 
explicitly as a diffusion equation by introducing the barycentric velocity as a reference velocity. The barycentric 
velocity v m is defined as 

v m = £>v, (2) 

Since the polymer mixture is assumed to be incompressible, we have that at every point 4>i — 1. (So far this does 
not exclude the possibility of a low molecular weight 'solvent' among the species i.) This mass conservation property 
implies that v m is solenoidal, i.e. V • v m = 0. If we introduce the following time derivative 

then the continuity equation, eq [l], can be rewritten as 

which has the form of a diffusion equation, expressed in terms of the diffusive flux of species i, which coincides 
with the total flux only when measured in the barycentric frame of reference: 

J i =J t i -p i v m = p i (v i -v m ) (5) 

We could, of course, have chosen any other reference velocity in place of v m , although in that case the continuity 
equation does not have the convenient form of eq ^. An exception is the special case of mechanical equilibrium 
(Vf> = 0); for this situation, the diffusive flux can be defined with respect to any reference velocity leading to the same 
expression for it. But in general, concentration gradients in our system will induce pressure gradients and thence fluid 
flow, which is handled by use of the barycentric frame. 

A. Phenomenological model 

Following Clarke's phenomenological model [Q, we can start from the diffusive (Cahn-Hilliard) 
dynamics expressed in eq ^ and assume that the diffusive flux of species i has the expression 

M° 

J, = -^V Ml + Pi v r (6) 
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where pi is the chemical potential of species i, and Mf is a phenomenological mobility coefficient, and v r is an 
unknown reference velocity that is common to all polymer species. We can use the 

incompressibility condition (local conservation of mass) to fix this reference velocity because J2i N-fi = 0. The 
diffusive flux then simplifies to 

Ji^^M.feV^fe (7) 

k 

where we identify the mobility matrix 

M° 

M lk = — J-fo - p t Ml (8) 

This model does not provide explicit expressions for the mobilities M?. ft seems reasonable to assume that they 
are proportional to the density, suggesting the more general form 

M ik = A t pi8ik - PiPkBik (9) 

where for length polydispersity At and Bik will depend on Ni. In general they might also depend on any other 
polydisperse feature such as the fraction of monomers of given chemical type on chain of species i. In the next section 
we develop a particular microscopic model in which the only dependences enter through the Ni (so that the case 
of chemical polydispersity becomes especially simple). We will show that we can recover the generic form for the 
diffusive fluxes, eq m, although now with explicit expressions for the mobility coefficients. 



III. MICROSCOPIC MODEL (FOLLOWING BROCHARD) 

For a dense mixture of long chains, we will use the tube model to describe polymer dynamics, as applied to 
interdiffusion by Brochard j|) . We assume that all the chains are dynamically similar in all respects apart from their 
lengths: that is, we ignore the effect of chemical species on entanglement length, mobility etc.. According to classical 
reptation theory (which assumes that each polymer moves through a tube made of the network) ||, the center of 
mass velocity v; of polymer species i, is related to its curvilinear velocity Wi (which determines its motion through 
the tube), as follows: 

Vi=Wi^- (10) 

where Li is the contour length of the tube, and is the end-to-end polymer vector. Under 
near-equilibrium conditions, these are related by 

U - ("J 

where the entanglement parameter N e is related to the tube diameter M. Note that N e should be smaller than N 
for all species present; for length polydispersity, this restricts the applicability of the present model to the case where 
all species, including any 'solvent', are fully entangled. 

In using eq [ll] we assume that the polymer conformation is not significantly distorted from its equilibrium counter- 
part. Hence, the dynamics we will describe below will hold only for chemical potential gradients weak enough that 
one does not enter the non-Newtonian flow regime, where the flow response to pressure gradients becomes nonlinear. 
Using the previous expression, the curvilinear and center-of-mass velocities measured relative to a static network are 
then related through 

\ v .\* = (12) 
N K ' 

In practice, the network of entanglements is common to all the polymers, but need not be static. If the system moves 
under the action of an external field, this network will be characterized by a common velocity, v t , or 'tube velocity'. 
Equation [l2| should then be generalized to 

|v,-v t | a = ^ (13) 

Following Doi and Onuki Jl0| , we may now deduce the (linearized) equations of motion by minimizing the energy 
dissipation rate of the system. This has two contributions: one related to the free energy variation, the second to the 
energy dissipation due to the friction between the polymers and the tube. We analyze the two terms in turn. 
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A. Free energy rate contribution 

Generically, the free energy of a polydisperse polymer mixture can be expressed as the sum of an ideal and an 
excess contribution: 

F = J drf[{ Pi (r)}} = J dr|fc fl T^ W (r)bgp i (r)+/ e!e |> i (r)]| (14) 

where f ex is in general a non-local function of the chain-number density. The free energy rate (which we define as 
the contribution to the dissipation of the system due to work done against thermodynamic gradients) can then be 
written as 

t£ = / dr = I ^ (r)v • {ViPi) dr (15) 

where the chemical potential is defined in general as /ij(r) = SF/Spi(r) and where use has been made of mass 
conservation, eq |l|. This result is general, although in section IV we will discuss specific examples for the free energy 
of polymer mixtures, expressed in terms of moments of the chain number densities. 

B. Tube dissipation rate 

The tube dissipation rate arises from the friction of the polymer against the network during its motion. Hence the 
relevant polymer velocity is the curvilinear one, wt. The dissipation is then assumed to be proportional to its square, 
and to the number of monomers of each polymer species. If we denote by £oi the microscopic friction constant of 
species i, we can write the tube dissipation rate as 

W = i drj^l^w* = f drJ2 UiWi - v t | 2 (16) 

where we have used the equilibrium relation between the curvilinear and the center of mass velocity for polymer 
species i, eq [l3[ and introduced the polymer friction coefficient, 

N- N 2 

6 = 6,^ = 6,^ (17) 

Since £ i is a microscopic friction coefficient, related to the interaction between a given monomer and the network, 
we assume it is the same for all the species, i.e. = While true for pure length polydispersity, when chemical 
polydispersity is present this will not be strictly correct. (The resulting errors should, however, be quantitative not 
qualitative in nature, especially, for example, when the chemical differences involve only deuteration.) 

Brochard's picture has allowed us to obtain explicit expression for the friction coefficients, derived from a microscopic 
model. In particular, it predicts a quadratic dependence of the friction coefficient on the polymerization index iVj of 
species i. Within our assumption of constant £oij it is only for length polydispersity that the mobility contributions 
will depend on polydispersity. For other situations such as chemical polydispersity, its influence on the diffusion 
coefficients will enter only through the thermodynamics of the system. 

Following Brochard [|| the tube velocity is now determined by requiring that the friction force acting on the network 
should balance. Such a force can be obtained from the dissipation rate W taking its derivative with respect to the 
tube velocity. Requiring 8W/ 8w t = determines the tube velocity 

v t = (18) 

Note that, except for the important case of length polydispersity, the tube velocity always coincides with the barycen- 
tric velocity in our model. This coincidence is very useful to analyse the diffusion, because under these circumstances 
Pi( v i — v *) is already the diffusive flux of species i (and we obtain below an explicit expression for Vj — Vt). In this 
respect, length polydispersity plays a distinctive role from the dynamical point of view. 
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C. Total dissipation rate 

The total dissipation rate (essentially the entropy production rate) for our polymer mixture is now written 

'1 



as 



R = fdr £ 



t&I v * ~ v t| 2 ~ MiV • (Pi v i 



(19) 



where the third term is introduced to enforce incompressibility (cf. eq ||). Here p, the associated Lagrange multiplier, 
corresponds to the pressure. 

The polymer velocities are derived by minimizing R. Imposing 5R/6~Vi = we get 



&(vj - Vf) + /OiV/ij + faVp = 
Summing over all components (and imposing incompressibility, i.e. X)i^( r ) = -0 we recover 

Vp = - y^p a V/Xj 



(20) 



(2i; 



which is simply the Gibbs-Duhem equation (this shows that the Lagrange multiplier p corresponds indeed to the 
thermodynamic pressure). Using this explicit expression, the velocity of polymer i has the form 



Pi <Pi \ 

v, = v t - — Vfii + — PkVp-k 
& 5i k 



(22) 



From this we can write the diffusive flux of species i as 



(23) 



We can again use the incompressibility condition to obtain an explicit expression for v t — v m ; imposing ^ Ni Jj = 
leads in this case to 



which shows that these two velocities are not independent of each other, and that their difference is proportional 
to the chemical potential gradients, but not simply to the pressure gradient (hence this difference is not caused by a 
departure from mechanical equilibrium). 

As we have already mentioned, in the absence of length polydispersity, v t — v m is zero within our model, because 
we assume £j cx 4>iNi. In the remaining case of length polydispersity, the difference can be expressed as 

v.-v^lW^-pW, (25) 
where p = p% is the overall local chain density, and a = £o/N e . 



D. Result for the mobility matrix 

Substituting eq |24] into the expression for the diffusive flux, eq [23], we get 

Pi c ( 4>i , 4>k 



J > = E I -f ** + '* (| + f J - p" E f | ^ ( 26 ) 

This diffusive flux is proportional to the chemical potential gradients. It has the form of eq J?], with the generalized 
mobility matrix, eq |9|. But in this case we have obtained an explicit expression for the mobility coefficients. As 
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had been hypothesized, the mobility is proportional to the chain number density. However, for this model we have 
derived the dependence on the polymer length (which enters through the friction coefficients and which turns out to 
be nontrivial). We can rewrite eq ^ in the same form as eq ^|, with the following choice of the mobility matrix: 

> f PiPi c j (<t>i .<t>k\ S^tf (V7\ 

M ik = — — 5 ik + pip k — + — - pip k > -7 (27) 

« V Si Sfc / j 

This shows that the phenomcnological model recovers the correct functional dependence for the diffusive flux. 

Comparing with the expression introduced in eq [| our analysis based on the physical hypotheses of Brochard's 
picture ||, yields 

A = - P j , S = -|i + |i + £f (28) 
Si Sfc . Sj 

In the subsequent sections we analyze in detail several applications of this result, for different forms of the polydisperse 
polymer free energy. 

As mentioned previously, our expression for the diffusive flux, eq has the same general form as that obtained 
by Clarke Q for the case of length polydispersity. The main conceptual difference is that, while Clarke expresses 
diffusive fluxes as gradients of the polymer chemical potential ( relative to the solvent chemical potential) plus a global 
convective contribution with respect to an unspecified common velocity, in our case we define diffusive fluxes with 
respect to the barycentric velocity. (We will come back to this comparison in section |vj.) So far, in fact, we have 
not introduced any explicit solvent, although the model as presented above does permit one to enter it as one of the 
species i, so long as this 'solvent' is itself made of chains that are long enough to be entangled. We return to the issue 
of solvent dynamics in the following sections. 



IV. FREE ENERGY FUNCTIONAL 



We now discuss the common structure of the free energy functional for all the systems we address below. For 
definiteness we consider polymer free energies that are derived from a specific nonlocal free energy functional proposed 
in ref 2. It has the form (setting fc^T = I) 



T = \ ^|^p fe (r)(log( Pfe (r))-f) + F-[{m„(r)}]| 



(29) 



where the excess free energy depends on a finite set of moments. These moments are locally varying quantities that 
depend nonlocally on the associated chain number densities through 

m„(r) = ^2 j N k w k (r - r')p k (r') dr (30) 

k ^ 

where Wk{r — r') is related to the structure factor of a polymer chain of species k in a specific way H]. The quantity 
Lo nk defines the nth moment in terms of the polydisperse variable; in all cases of interest here it is a simple power 
(e.g. uj n k = in the case of length polydispersity). 

The use of this functional, which contains the information on the internal structure of the chains, will allow us 
to address issues such as the dependence of the initial unstable modes on the wave vector, in the case of a system 
undergoing spinodal decomposition. For small wavevectors the nonlocal kernel ujk appropriate to a Gaussian chain is 
well approximated by the following gradient expansion 

w k {r - r') = (1 + ^V 2 )<5(r - r') = W Qk S(v - r') (31) 

where the additional subscript denotes the small wavevector limit. The operator Wqu will be used in subsequent 
sections. Note that the chosen normalisation, J u>fc(r) dr = I, means that factors N k appears beside it in eq |3C|. 
Thus, in the case of length polydispersity, the zeroth nonlocal moment mo(r) (which corresponds to choosing wofc = 
N® = 1) in fact describes the local monomer concentration. For chemical polydispersity it is again the local monomer 
concentration, with no discrimination between monomers of different chemical species. (In the absence of a solvent 
species mp(r) is unity everywhere, because of incompressibility.) 
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Note that although we use a specific form, eq £9|, for the nonlocal contributions to the polymer free energy, it is 
chosen so as to reduce to a standard Flory Huggins theory, in the limit of zero wavenumber. This requirement fixes 



the form of the term F ex [{m„(r)}] in eq 29 Moreover, once the gradient expansion eq 31 is adopted, the nonlocality 
we choose has a wider justification than the underlying model defined by eq |3^. In particular, the same form for 
the gradients could be obtained by the random phase approximation H]. A non-local free energy functional will in 
principle require, for consistency, a non-local expression for the mobility matrix. Such effects have been considered 
previously in the literature p^ |. However, for simplicity, we will make in this paper the assumption that the mobility 
coefficients can be regarded as local, as is normal practice || f|, 

V. CHEMICAL POLYDISPERSITY 

Chemical polydispersity describes a set of polymers all with the same length, but a continuously varying interaction 
parameter. The normal example is a melt of random copolymers of equal length made of two monomers, e.g. 
deuterated and hydrogenated, that differ slightly in their enthalpic interaction. Within the model outlined above, we 
assume that this is the only difference. Each chain is then described by a polydisperse quantity fi obeying — 1 < / < 1 
which controls the proportion of the different monomers present along it; they have no correlation other than this, 
and within a mean field approximation the presence of variable / is equivalent to a spectrum of x parameters affecting 
different chains |Ts) ]. The moment densities m n (r) are then defined by the choice aj n k = f^\ and within Flory Huggins 
theory, the excess free energy is simply F ex — —xm\ [fjj. Here \ is the Flory parameter controlling the interaction 
between the two different monomers. 

Next we treat this case of the chemically polydisperse melt. However, we also allow for a solvent to be present. 
Technically, this is already included in the preceding formalism so long as the solvent is considered as 'just another 
species' within our incompressible mixture of many polymers. Indeed there would be nothing to prevent us developing 
the general formalism to address this case without singling the solvent out for any special treatment at all. But the 
fact that its length and chemistry will usually be different from all the remaining polymers implies that the solvent 
will have a distinct dynamical behavior, which it is convenient to handle separately. Moreover, because of 
incompressibility, both the thermodynamics and the dynamics of the solvent are subservient to the remaining species 
and this will allow us to focus on the dynamics of the moments m„(r) of the chemically polydisperse chains alone. 
These dynamics will, of course, be different according to whether a solvent is present or not. 

Note that this treatment of the solvent does not release us from the assumptions of the tube model made in our 
earlier derivation of the diffusive fluxes. In particular, this means that the solvent itself must consist of chains long 
enough to be fully entangled (though, in the cases addressed below, much shorter than the remaining, chemically 



polydisperse, chains of primary interest). However, as we will discuss in section VI, we can later on use the general 



form of the equations for the diffusive fluxes to guess how these may behave for an unentangled solution. 



As discussed in section [II above, a key simplifying feature for chemical polydispersity is that (within the approxi- 
mations we have taken) the structure factors and mobilities of the chains are independent of the polydisperse feature 
fi. In fact, many of the results described below for this case will also hold for any other kind of polydispersity with 
the same attribute, even if the structure of the excess free energy is not that corresponding to chemical polydispersity. 
We will point out, below, where the number of moments appearing in the excess free energy plays a role in the results 
we obtain. 

A. Chemically polydisperse polymer melt 

The diffusive flux for this system is specially simple. The mobility coefficients in eq |2^ reduce to 

A = ~ ( 32 ) 
aN£ 

b = -^r ( 33 ) 

With our assumptions, this is true for any melt when all polymers have the same length N = N p . 

We will focus on polymer melts whose excess free energy depends only on a finite set of moments. With chemical 
polydispersity, Flory-Huggins theory states that only the first moment of the chemical composition is involved ED]: 

F ex { mn (r)} = -xml (34) 
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Correspondingly, the chemical potential gradients simplify to 

\7p k = — \7p k - 2xN p u lk 
Pk 



1 + — ^V 2 
24 



Vtoi 



(35) 



Note that in this case the function Wo = 1 + -gjV 2 does not depend on the polymer species. Since all polymers have 
the same length, we have (to second order in gradients), 

N p ^u nk p k (r) = 



N v 9 

1 ^V 2 

24 



m„(r) = W 1 m n (r) 



(36) 



where W Q 1 is the inverse operator of Wq, likewise expanded to second order in gradients. Making use of this relation, 
the evolution equation for the species reads 



dpi L-Y^v-f-— 

dt ~ aN n 4^ V N„ 



Pi PkS P-k 



(37) 



where, due to incompressibility, Vmo = 0. This equation holds to leading nontrivial order in separate gradient 
expansions in the thermodynamics and in the relationships between fluxes and chemical potentials [12]. This is the 
normal procedure in Cahn-Hilliard approaches to polymer dynamics and is the simplest theory able to predict a 
nontrivial wavevector dependence for the growth rates in spinodal decomposition |5| ^, [To| |. 

The evolution equation of the moment densities TOj(r) can in turn be written, within the same order of approxima- 
tion, as 



dm r , 
~dT 



= D Q {V 2 m„ + 2 X N P W Q V ■ ((W^mi)^ V0 - {W^m^WoVrm} 



(38) 



with Do = l/(aN 2 ). The diffusion coefficient Do sets the time scale, which as expected, decreases proportionally to 
the increase of the friction coefficient. 

This equation for the moment dynamics can be easily generalized to the case in which the excess free energy is a 
function of any finite set of moments. One gets in that case 



dm n 
dt 



= Do 



V 2 m r - 



N P W V -J2 F af ((Wo- 1 ™^) - {W^m^W^m^) W o Vm 

a(3 



(39) 



Note that for generic free energies this expression involves a linearization in moment gradients (which was not required 
for an excess free energy quadratic in mi as considered above) . Apart from this, eqs [38] and [39] have the same structure. 
Each moment relaxes proportional to its own density deviation with a relaxation rate Do, but is also coupled to all 
other species through the various moments that appear in the excess free energy of the model. 

Note also that the mobility does not play a major role; it only defines a characteristic time scale (which we will call 
the fast scale) Tf(q) = 1/D q 2 where q is the wavenumber. In contrast, the coefficients F£% are related to the location 
of any phase transitions and/or spinodals. Proximity to these will control the form of any slow modes present. As a 
result, for an excess free energy with n moments there are n potential slow modes related to the thermodynamic forces 
of the system. This implies that, even if there is an hierarchy of kinetic timescales, it will be finite. (The situation will 
be quite different for length polydispersity, where there is no clear factorization of the moment evolution equations 
— see section VI below.) 

For the chemically polydisperse melt there is only one moment appearing in the excess free energy (mi). The 
eigenf unctions of the diffusion matrix|^l| then have a very simple structure. In fact, all moments with n > 1 are 
eigenvectors with the same eigenvalue A„ = —q 2 D |3 with no wave vector dependence other than the usual diffusive 
one. On the other hand, when we linearize the evolution equation ^ at wavenumber q we find the one remaining 
eigenvalue 



Ai = -q 2 D [1 - 2 X N p (m 2 - mfW^q)] = -D Q q 2 



X 



l--(l-'-^q 2 ) 
Xs 12 



(40) 



where the tildes denote values in the homogeneous parent state about which the linearization is performed. We have 
also used that the interaction strength at the spinodal is given by N p \ s = l/(2(m2 — fh 2 ))- As expected, Ai changes 
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sign at the spinodal, and defines a slow mode in its vicinity (slow relative to all higher eigenmodes A„>i) that will 
control the early-stage kinetics of the phase separation process. During the initial stages of phase separation all 
moments n > 1 will relax fast toward their corresponding local equilibrium values; but these values are controlled by 
Ai which evolves on a slower time scale. Note that this slow mode is not simply the first moment mi, but a linear 
combination of all moments. 

These results for the chemically polydisperse melt contrast strongly with those given in the previous literature for 
length polydispersity [j| as discussed further in section VI below. 

Specifically, since there is only one relevant moment, there is no distinction between a quenched and an annealed 
spinodal that can define different kinetic scenarios as described by Warren jj| for length polydisperse solutions. More- 
over, the slowest mode (which is the unstable one) is the only mode that involves the first moment. All moments of the 
distribution, except for the first one, will accommodate rapidly to the mean local composition, governed by the slow 
relaxation mode involving mi . Finally, the wave vector dependence of the unstable mode is independent of polydis- 
persity, and the same as in a monodisperse system. Only the overall amplitude (determined by the thermodynamics) 
is modified by polydispersity. 



B. Chemically polydisperse polymers in a passive solvent 



In this case, unlike the melt, the mobility plays a role beyond that of a simple scale factor in the dynamics. As 
explained previously, it is convenient to single out the solvent for special treatment en route to deriving dynamical 
equations for the moments m n . For dcfinitcness (and a clearer separation of timescales) we treat the case where the 
solvent has length N s <c N p , although we continue to assume that all species are entangled. The ratio of the two 
lengths is defined as r = N p /N s . 

The generic expression for the diffusive flux and the free energy will both depend now on the solvent number density. 
We first rewrite the relevant expressions in terms of polymer densities alone. These expressions will be valid both 
for this model and for length polydispersity addressed in section VI. Due to incompressibility, the solvent number 
density is related to the monomer concentration (mo) through 



Pa 



1 - mo 



(41) 



where we have disregarded the internal structure of the solvent in writing this as a local relation. This is a good 
approximation at large r when the solvent is much smaller than the remaining polymers. In what follows, we treat 
r as a large parameter and are mainly concerned with the leading behaviour in that limit; but where subleading 
corrections are calculated, this is done without correcting eq EI] , 

We can now single out the solvent contribution both to the diffusive fluxes and the chemical potential gradients. 
For what we will call a 'passive' solvent, its contribution to the free energy, eq is simply to add a local entropic 
term J o s (lnp s — 1). Then 



Ps^Ps 



1 



■ Vmo 



(42) 



Note that for symmetric distributions in a passive solvent there is no tendency of the chemically polydisperse polymers 
to separate collectively from the solvent; the only tendency towards phase separation is within the polydisperse 
distribution itself. This ceases to apply for an asymmetric distribution of solute species. Note also that a more 
general solvent would introduce (as well as ignorable linear terms) terms in momi into the excess free energy, leading 
to somewhat more complicated phase behavior Q and dynamics. 

Using the fact that p s depends only on the overall monomer concentration, the diffusive flux for the polymer species, 
given by eq ^, can now be expressed as a function of polymer variables only: 



4>) (1 - m ) 5 



— Vl7in 



Pkpi^pk 



(43) 



where the ' means that the sums run only over the solute polymer species. 
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For a solution of chemically polydisperse polymers, the excess free energy for a passive solvent may again be written 
within Flory Huggins theory as 

F ex = - X m\ (44) 

where the interaction amplitude of interest (around the spinodal) obeys X N p — 1. The chemical potential gradient is 
again given by eq |3^. We can now derive expressions for the diffusive fluxes of each polymer species, following the 
same kind of derivation as for the melt. We take advantage of the fact that the friction coefficients for all species is 
the same (see eq |l7]) to arrive at 



J 4 = D pi ^2 {Sik + N p p k } V/Xfe 
k 

- D oA |iV p (l-mo)(r-l)^p fe VMfe + r(l + mo(r-l))Vmo| (45) 

The first line coincides with the expression of the diffusive flux of a species in the melt. The presence of a passive 
solvent has two effects in the polymer dynamics, which correspond to the two terms appearing in the second line of 
eq ^5]. Gradients in the local monomer concentration generate a flux of each species of polymers, but in addition, the 
local presence of a solvent allows global polymer rearrangements according to 

' 1 

$>fcV Mfe = — W^Vm - 2 X {W^ 1 m 1 )W Q Vm l (46) 
k ly P 

The analysis is very different from the melt case because now gradients in the local overall monomer concentration 
do not vanish. 

Introducing the diffusive fluxes of the moment densities 

j„ = N p J2^nkW Q J k (47) 

k 

then using eq ^5] these read 

j„ = D {-Vm„ + 2 X N p W [(W^m^) - (W ~ 1 m„)(W - 1 m 1 )] W V mi 
+ 2 X N p W (l - m ){r - l)(W^ 1 m n )(Wo 1 m 1 )WaVmi 

+ Wo{Wo 1 m n )[^r{l-m +m r) + (2-mo-r(l-m ))W ' 1 ] Vm } (48) 



leading to 



D {-Vm„ + 2 x N p [m n+ x - m„mi] W 2 Vmi 



+ 2 X N p (l - rh Q )(r - l)m„TO 1 W Vm 



+ rh n [2 — too — r(l — mo) — r(l — too + TOor)Wo] Vtoo} (49) 

where we have now linearized the fluxes, as appropriate for studying the structure factor during the early stages of 
phase separation. As in eq EoL the tilde refers to values of the moments evaluated at the reference state, instead of 
being dynamical quantities. 

Again, the first line in the previous expression is the contribution that we obtained for the melt. The second line 
requires the presence of a solvent and corresponds to a correlated flow of all the species that avoids segregation. 
Finally, the third line corresponds to the diffusive flux of a given species induced by solvent gradients. 

In terms of these moment fluxes, the evolution equation for the moments read 

drrin V7 • rxris 
~dT = in ( ' 

which yields the final results for the linearized dynamic equations of the moment densities themselves 
^ = Z? {V 2 TO„-2 X iV P [TO„ +1 -TO„TO 1 (l-(r-l)(l-TOo))]PFo 2 V 2 m 1 

- m„ [1 - (1 - mo)(r " 1) " Kl + rh (r - 1))W ] V 2 to } (51) 

Although these equations are more complicated than for the case of a melt (eq Q), they do retain a relatively 
simple quasi-triangular form |l5| which allows us to analyze in some detail the structure of the eigenvalues. The fact 
that all moments only couple to the first two moments (too and mi) imply that, except for these two, the rest are 
degenerate eigenfunctions. Hence we may set A„ = —D$q 2 for n = 2, 3, ... . The interesting dynamical behavior lies 
then in the remaining two eigenmodes, which are the only ones that involve toq and mi. 



11 



1. Dependence on solvent size ratio 



In general the expressions for the two eigenvalues is quite involved, but given our assumptions so far, they can be 
found for arbitrary r. However, the derivation was subject to eq ^ which itself is valid for large r only. Hence we 
relegate the general expressions to Appendix Hi 

In the limiting case of very large polymers (i.e. r»l) the two eigenvalues simplify to 



Ao 
Ai 



-q 2 D r 2 rh 2 W (q) 



q 2 D 



1 



Xs 



-W Q (q) 2 



(52) 
(53) 



where Xs = mo/ (2(mo?ri2 — to?)) is the spinodal value of the interaction parameter (to leading order in large r). For the 
particular case of a symmetric parent (i.e., a symmetric distribution of the polydisperse species in the homogeneous 
phase prior to the quench) the expressions for the the eigenvalues are particularly simple. In this case fhi — in the 
initial state and the eigenmodes that involve toq and mi do not mix them; their eigenvalues (exact to order q ) are 



Ao = -D q 2 [1 - m {1 - (1 - mo)(r - 1) - r[l + m (r - l)]W (g)}] 

.2 (■> n ,~ itt f„ \2\ r> _2 / 1 X 



Ai = -D Q q* (1 - 2 X m 2 Wo(qy 



-D Q q 2 1 - ^W (q) 2 

Xs 



(54) 
(55) 



with Xs = 1/(2to2) the spinodal value of the interaction parameter for a symmetric parent. Note that eqs ^ and 
|55| have the same form as the relevant eigenvalue in the melt case, eq The only difference is the location of the 
spinodal; it is through the thermodynamic sector only that the solvent, and likewise the polydispersity, can enter. 
This shows that in the limit of high size ratio (r 3> 1) the main features of mode structure are, within our model, 
universal. 

Likewise if we look at the wave vector dependence of the eigenvalues Ao and Ai, we see that it does not depend 
directly on polydispersity. This behavior is again analogous to that of the melt, and can be attributed to the absence of 
coupling of the unstable mode to all the rest of the modes. This is contrary to what happens for length polydispersity, 
as we discuss in section jv| . 

For both symmetric and asymmetric parents, in the limit r»l, the two eigenvalues scale asymptotically as 



Ao 
Ai 



2-2 

q m 2 

'aN 2 

aN 2 



1-2L 

Xs 



(56) 
(57) 



but the convergence to this limit as r is increased strongly depends on parental asymmetry. In Figurel we show the 
wave vector dependence of the unstable eigenmode, Ai for an exponential parent (p n ~ exp(a/„)) as a function of 
the size ratio r for different values of a. When the asymmetry parameter a is increased, the wavevector dependence 
becomes more sensitive to the size ratio r (in fact, for a — there is no dependence on r; cf. eq [35|). The dependence on 
size ratio becomes also more relevant as the polymer dilution increases. For high asymmetry and dilution (Figurelc), 
the convergence to infinite aspect ratio becomes much slower and it is not monotonous for intermediate values of r. 
In Figure2 we show the wave vector dependence of the stable mode, Ao, for two disparate values of the asymmetry 
at high dilution. The dependence on q is monotonous as anticipated. In all cases the curves converge in the limit of 
large aspect ratios (note that because of the different scaling of the two eigenvalues Ao and Ai we have scaled here 
Ao by r 2 ). For large asymmetry (Figure2b), also for this stable mode the convergence to the high aspect ratio limit 
becomes slower and it is non-monotonous for intermediate r values. 



2. Qualitative behaviour at large size ratio 



The above analysis shows that, for large size ratios between solute and solvent (r ;§> 1), Ao is much larger than the 
rest of the eigenvalues. This means that the mode connected with Ao (which, for a symmetric parent only, does not 
involve mi) will grow or decay much faster than all the rest of the modes, whose rates scale with the inverse of the 
solute polymer length N p , defining for n>la common time scale Tf(q) = 1/D^q 2 . The latter was already identified 
as the fast time scale in the melt case. The corresponding modes remain fast relative not to the solvent mode (which 
is faster), but to a single slow mode. This is governed by Ai, which splits off from this family of modes with a time 
scale that diverges on the spinodal. 



12 



This means that during the early stages of phase separation (or during the decay of an external perturbation 
within the stable region of the phase diagram), the overall monomer density is first equilibrated (if necessary); on 
this timescale the remaining moments are effectively quenched. In a second stage the chemical composition will relax 
towards its equilibrium distribution (with the monomer density slaved to this). Close to the spinodal this chemical 
relaxation can, just as in the melt case, be separated into two stages. These entail relatively rapid equilibration of 
the higher moments to a state of local equilibrium set by the slow evolution of the mean composition variable mi 
(which is the local excess of A over B monomers in the case of a random copolymer system). It is notable that, 
for an asymmetric parent, the two relevant modes include nontrivial linear combination of mo and mi, while for a 
symmetric parent these two moments appear separately in the corresponding eigenfunctions. This is because of the 
thermodynamic coupling which, for an asymmetric parent, requires chemical separation to be accompanied by changes 
in the overall monomer density. 

This picture is somewhat analogous to one described by Warren for the case of length polydisperse polymers 
undergoing phase separation from an incompatible solvent Warren developed a picture in which there were two 
separate spinodal curves: the equilibrium one and, within it, a quenched one to describe the dynamics of a hypothetical 
system whose overall monomer density could change but where the relative prevalence of different chain lengths could 
not change locally. Between these two spinodals, the kinetics were argued to be controlled by the slow sorting out of 
the polydisperse species. 

For chemically polydisperse chains in a solvent, a similar scenario can arise. However, for the passive solvent 
considered here, this can only happen for an asymmetric parent since otherwise the phase separation is towards two 
states of equal solvent density; this density coincides with that in the initial state. Accordingly there is no driving 
force for phase separation in a symmetric system on the time scale of the fast (solvent) mode as may be confirmed by 
constructing the relevant free energy with quenched (/) = rhi/rho. If we consider a fully quenched situation where all 
higher moments are slave to mo, ( i.e. higher moments evolve as m n = (/™)mo), its linearized evolution equation is 

^ = A>{l-2 X mo</> 2 (l-mo)[l + mo(r-l)]^o 2 

+ ?no[-l + (l-m )(r-l) + r(l + (r-l)mo)Wo]}V 2 m (58) 

which for large r will give a relaxation rate scaling with the fastest mode of the system. This leads to a rapid growth 
of fluctuations, but only inside a region where the quenched system is unstable; such a region is present only for an 
asymmetric parent. 

In Figure3 we show the quenched and annealed spinodals for a chemically polydisperse mixture with an exponential 
parent shape (p n ~ exp(a/„)) for different values of the parameter a (which controls the degree of asymmetry of the 
mixture) and of the size ratio r. One can see that in the regime where polymers are much larger than the solvent, the 
two spinodals are far apart. There is thus a broad region, in between the two spinodals, where the density relaxation 
will take place only through a slow sorting of the chemical species. The slowest process will then control the initial 
steps of the phase separation process, determining the relevant time scale. 



VI. LENGTH POLYDISPERSITY 



A. Structure of equations 

Before analysing the effect of length polydispersity in detail, it is instructive to compare the structure of eq ^6| with 
the analogous one derived by Clarke M. He obtained diffusive fluxes for the monomer concentrations, disregarding 
(as we do) the internal structure of the chains when constructing the relation between fluxes and chemical potential 
gradients. (The internal structure enters only when calculating these gradients.) We may rewrite our results in terms 
of the monomeric flux of each species i (that we will call Jj) starting from eq ^6|. We get, 

J* = E j -f + (| + |) - E f J tiWfiu - E A ^ ( 59 ) 

where jlk is the monomeric chemical potential of species k 7 i.e. jlk = SF/S<f>k- The mobility matrix A^ coincides with 
the one derived by Clarke (cf. eq 22 of Ref. Q) if we identified the phenomenological coefficients Xj (l(| introduced 
in that paper as 



1 6 



(60) 



13 



The basic differences between Clarke's approach and ours is then that while he considered diffusive fluxes referred 
to a common velocity fixed through incompressibility, in our case the diffusive fluxes are defined with respect to 
the barycentric velocity. In this way we get diffusive fluxes that are strictly proportional to the chemical potential 
gradients of each species. 

Nonetheless this comparison is suggestive, because the phenomenological model of Clarke is claimed to apply both 
to entangled and unentangled mixtures. In the former case he suggests that Ai should scale as the inverse of the 
polymer length, A.; ~ 1/iVj. We have derived this result in section III. Brochard's theory (which was our own starting 
point in section II] ), and also Kramer's theory [|l7) (on which Clarke based some of his reasoning) both assume that 
the diffusion is controlled by the fastest species (they are accordingly referred to as 'fast mode theories') |l8| |. Hence, 
it is not surprising that we get mobility coefficients that concur with those posited by Clarke on the basis of Kramer's 
model. 

For the unentangled case, Clarke argues that, because of the different scaling of the self-diffusion coefficient with 
polymerization (Rouse model ||) one can assume that in this situation A^ is a constant instead. Although we have 
derived the expressions for the diffusive fluxes within the context of the tube model, if we assume, following Clarke, 
that their overall functional form is the same in the unentangled case we can use our results in both situations 
simply by changing the dependence of the friction coefficients on the polymer lengths. For chemical polydispersity as 
discussed above, this means modifying only one or two global parameters (solvent and s olute chain mobilities) of the 
theory; for length polydispersity the effect is more complicated and discussed in section VI D| below. 



B. Length-polydisperse chains in a solvent 



For length-polydisperse chains in a solvent, treated within Flory Huggins theory, the free energy is written as 

T= Adr5^pk(r)(lnp fc (r)-l)+xmo(r)(l-mo(r)) (61) 
J k 

where we have initially chosen the sum over species to include the solvent chains, while the moments are the linear 
nonlocal combinations of the polymer number densities defined in eq [30] with u n k — NJ}, and the solvent chains 
excluded. Note that, as described following eq|3^, mo is the total volume fraction occupied by solute monomers. The 
chemical potentials of the species are accordingly 



/i fe = log p k + xNk - 2xN k 



, J w k {r - r')rao(r') dr' 



fJ-s 



\ogp s = log 



1 - rap 



(62) 
(63) 



where we again neglect the spatial extent of the solvent chains, presuming them small compared to the typical size of 
the solute polymers. 

As we did for chemical polydispersity, we now choose to treat the solvent as a separate species and derive equations 
for the diffusive fluxes of the solute polymers with the solvent coupling eliminated (statically and dynamically) via 
the incompressibility constraint, eq ^J. Using this relation, and the fact that Woi — 1 is linear in species length p9L 
the diffusive flux of species i can now be found to second order in gradients. The full result is given in Appendix^ 
and is quite complicated. Truncating further to first order in gradients, the fluxes obey (with a = £,o/N e ) 



1 

N. 



2 Vpi + piVm_ 2 



+ Pi 



2 X (l- mo )-- 



Pi_ 
1 



TO_i 



1 - mo 



too . 

iv7 |Vm ° 



(64) 



The completely different appearance of this expression from the one obtained previously for 
chemical polydispersity arises because the friction factor £ in eq |26| depends explicitly on 

polymerization index Ni which is now the polydisperse variable and not the same for all solute chains. Because of the 
factor Nk arising in the definition of the moments, eqpfl, the overall number density of solute chains is m_i. Since 
in addition the friction appears as a denominator in eqpq, two negative moments appear in the above equation. 
Rewriting this in terms of diffusive fluxes j n for the moment densities m n (r), we obtain in the linearized regime 



ajn = -Vto„_ 2 + to„Vto_ 2 



m n -i - m n m^i - m n - 



1 - to 



N, 



Vto_i 
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2x(l-mo)-^- 



m„_i - m„m_i + m °™ n ) Vm i ! 



In Appendix |b| we provide the corresponding expression to second order in gradients, which can be used, in principle, 
to study the wavevector dependence of the relaxation modes. 

But even in the long wavelength limit, the set of equations is no longer at all tractable: moment m„ is coupled to 
m n ^2, and all moments are coupled to too, to_i and m~2- The higher moments are no longer simple eigenfunctions. 
Also, there is the additional complication that moments of negative powers of the chain length have entered the 
description. Since we know that the set of functions u>i(N) = N l form a basis for i = 0, 1, we can in fact express 
each negative moment as a linear combination involving only moments with positive index. In doing so, though, we 
reach an expression for the diffusive fluxes of positive moments which couple to gradients of all the other positive 
moments. It is clear that the diffusion matrix is no longer truncatable and an infinite hierarchy of dynamical modes 
can be expected. Note that this hierarchy remains compatible with a gross separation into two basic timescales when 
solvent and solute chains are very different in length. One can still expect a single mode scaling with the inverse 
solvent size, with the remainder involving the solute chain sizes instead; the hierarchy controls the details of how the 
solute chains segregate by size. 



C. Approximate truncation of the dynamics 



Following Warren j^] , we can try to analyze in more detail the structure of the dynamic modes for length polydisper- 
sity by considering the dynamical equations of a finite subset of moments created by simply truncating the diffusion 
matrix. Thus we assume that the flux of any species is only coupled to spatial gradients of a finite set of moments. A 
priori, there is no guarantee that such a process will preserve the character of the modes; we will attempt to verify, a 
posteriori, that this procedure is in fact sensible. 

Let us consider a two-moment description and assume that the two relevant modes are the volume fraction, too, and 
the chain number density, m_i. According to eq|6^, their fluxes are coupled to gradients of other moments. However, 
the assumption that only too and m_i are relevant allows us to express other moments by the corresponding projection 
into the subspace spanned by the retained moments. To be precise, the weight functions LOk with k different from 
and —1 will be approximated as a linear combination of those two weight functions, using the shape of the parent. In 
this way, we can write 



{N- k ){N- 2 } - (N-^iN- 1 ) (A- fe - J ) - (N~ k )(N~ 1 ) 



(N- 2 ) - (N- 1 ) 2 



-w - 



(N- 2 ) - (N- 1 )' 



(66) 



where (...) means an average over the chain length distribution in the parent. With this choice we ensure that the 
average of and w-fcA -1 over the parent remains correct. 

We can now approximate the two additional moments (to^2 and 777^3) appearing in the evolution equations of the 
monomer and chain number densities as linear combination of the retained quantities. In this way, we disregard the 
coupling of all moments except for too and m_i, getting a closed set of equations. The required projection of moments 
to_3 and m_2 read 



TO_ 3 



TO_ 2 



-toq 



(N- 3 ){N- 2 ) - (N-^jN - 1 )^ t (N- A ) - {N-^iN- 1 ) 

(N- 2 ) - (A- 1 ) 2 
(A- 2 ) 2 - (A- 3 )(A- 1 ) 



(A- 2 ) - (A- 1 ) 5 



too + 



(A- 2 ) - (A- 1 ) 2 
(A- 3 ) - (A- 2 ) (A- 1 ) 



m -i = fhl m -i + Pm m a 



(A- 2 ) - (A- 1 ) 2 



m -l = P2l m ~l + 02O m O 



(67) 



The averages (. . .}can be expressed in terms of the moments of the parent (which are constants, not dynamical 
variables) through (A n ) = to„_i/to_i. 

The diffusive fluxes for the retained moments then reduce to 
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(68) 
(69) 
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-/?2i + m_i - 



ajo = (1 - mo) 



+ (1 - m ) i -fa + 2 X 



too 
N, 



m_i(l - toq) + 



m_i 



(1 - m )iV2 



V?71o 



(70) 



In the linearized regime, each of the factors preceding the gradient operator can be evaluated using the parent, in 

/ Arfc-f-l\ 

which TOfc = toq jj^k . Using also the size ratio r = (N)/N s , the linearized equations read then 



aj-i 



-(N) 2 (3- 3l + (N)(3- 2l m + ^((AT 1 )^) - m - r(l - mo))] 



(TV) 
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(71) 



These expressions remain fairly cumbersome. To proceed further and gain explicit expressions for the eigenvalues, 
we will consider only one specific example, where the parent has a uniform distribution with mean length (N) and 
half width A. We will work with the preceding expressions which are lowest order in the gradient expansion. To 
simplify things further we will consider the limit of a narrow parent ({N) 3> A) and large size ratio (r 3> 1). To 
lowest order, the two eigenvalues then simplify to 
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(72) 



where D m = l/a(N) 2 is a characteristic diffusion coefficient for the polymer solution, and where the value of the 
spinodal strength \s coincides with its monodisperse counterpart to the order of validity of the previous expressions. 

We see that (at least in the low q limit studied here) the two modes show a structure analogous to, but crucially 
different from, that described before for chemical polydispersity. The generic tendency at large size ratio r is for Ao to 
greatly exceed Ai, with the former defining a fast time scale and the latter a slow one. But this tendency is reversed 
close enough to the spinodal, where Ao now vanishes (not Ai as in the chemically polydisperse case). The difference 
arises because for chemical polydispersity the excess free energy involves the first moment of the polymer species, while 
here it involves the monomer concentration (zeroth moment). The nominally 'fast' mode is the thermodynamically 
unstable one; this was anticipated in Warren's original proposal Q. Between the quenched and the annealed spinodals, 
the phase separation can only proceed initially via species sorting, whereas beyond the quenched spinodal, the fast 
mode (Ao) will drive the phase separation process from the outset. 

In order to see if the assumption that only two moments matter is too crude, we have repeated the analysis assuming 
that the relevant moments are the three to which all the other moments are coupled, i.e. to_2,to_i,too. For the 
linearized regime we give the expressions for the corresponding fluxes in Appendix O In the limit of a narrow parent 
and large asymmetry the qualitative behavior of the three eigenmodes, found numerically, is the same as the one 
discussed with the previous approximation , in the sense that in both cases one eigenvalue (Ao) is generically larger 
than the rest (since it scales with the inverse solvent size), but crosses over to become smaller than the others in a 
narrow region as the spinodal is reached ]2C||. 
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D. Unentangled solution of length-polydisperse chains 



As we have discussed previously, if we follow Clarke B , we can assume that the expressions for the diffusive fluxes 
derived for an entangled mixture in eq ^6] keep their general form for an unentangled solution. The only difference is 
that in this case the friction coefficients are proportional solely to the monomer concentrations, £j = ot<f>i, and beyond 
this have no explicit dependence on polymer length. In this case the diffusive fluxes of species, eq reduce to 



J, 



and those of the moment densities become 
1 
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where we have related local moments of chain number densities to local moments of monomer 
concentrations fL9| . 

To analyze the onset of phase separation, we can linearize the previous expression to give 



■Vm„_i + m„Vm_i + rh„ [2x(N)(l - m ) - r 
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If we assume that we can isolate the dynamics of mo and m_i by closing their dynamic equations as done in the 
previous section, we obtain for a parent with narrow uniform distribution, in the limit of a narrow distribution (i.e. 
A(l) 
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where, as before D m = l/a(N) is a characteristic diffusion coefficient for the polymer solution. Again, we have kept 
only the dominant contribution in r. The value of Xs is the same as for the entangled solution, and deviations from 
its monodisperse value are also negligible in the narrow limit examined. 

The eigenmodes here are thus essentially equivalent to the situation described above for the entangled solution. 
We again observe the appearance of a 'fast' mode which nonetheless becomes slow very close to (and then unstable 
at) the spinodal. Away from there, the generic ratio between the two eigenmodes in this case is not as big as for the 
entangled case, because the dependence of the friction coefficients on polymerization is not as strong here. 

Incidentally, note that in the hypothetical case that the friction coefficient would have been proportional to the 
chain number density, then all moments would have coupled to themselves and to mo, mi and m-2. In this case, the 
only nontrivial eigenmodes would have been the three involving these three moments. A simple picture in terms of 
the moments is hence recovered exactly in that case. 



VII. CONCLUSIONS 



In this paper we have analyzed the dynamics of polydisperse polymeric materials. Polydispersity affects both the 
thermodynamics and the mobilities of the system. For a polymeric dense solution, for any polydispersity other than 
length polydispersity, the mobilities are the same for all species. In such case the moment structure of the excess 
free energy carries over to the dynamics of the system, and it is possible to obtain a simple picture in terms of the 
moments that appear in the excess free energy. Those are the relevant ones in the linearized dynamic regime, and we 
have analyzed their implications for the particular case of chemical polydispersity. 

The moment structure of the equations of motion makes it possible to map the polydisperse system into an effective 
binary mixture with a judicious choice of "effective" species. These fictitious species are linear combinations of chain 
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densities such as the monomer concentration and the mean chemical composition. Previous work has not exploited 
this choice but instead used a small number of discrete values of the polydisperse variable to impersonate the case of 
continuous polydispersity 0. 

We have seen that this finite moment structure implies in turn that there exists a finite hierarchy of relevant modes 
in the system. These modes define a fast and a slow time scale that will affect how kinetics takes place depending on 
the position in the phase diagram. Hence, there will be situations where the fast separation of the overall density will 
be dominant at short times (and only later the local chemical composition will rearrange on a slower time scale) , but 
in certain instances it will be the chemical sorting of species which will control the kinetics from the beginning. 

The difficulties in the treatment of length polydispersity come from the interplay of dynamics and thermodynamics. 
In the approach of Warren || such a coupling is not accounted for. Were no such coupling present, our analysis could 
follow that for chemical polydispersity; we might even be tempted to average (over the parent) the diffusion matrix and 
the chemical potential gradients separately. However, our analysis suggests that this constitutes a poor description of 
the true dynamics. 

The dependence of mobility on polymer length makes the treatment of length polydispersity more complicated due 
to the interplay between free energy changes and mobility couplings. Starting from a fast mode theory, we have derived 
the equations of motions for the moments. By keeping both positive and negative powers of the length, we arrived at 
a set of dynamical equations that, while more complicated than for chemical polydispersity, retain a relatively simple 
structure. In fact, the evolution of a given moment is only coupled to gradients of three other moments (whereas 
in terms of the species, all species are coupled to gradients of all the other species). The fact that we have to keep 
negative powers of JVj stems from the fact that we are dealing with a fast mode theory. If we re-express those powers in 
terms of positive powers of JVj, then the dynamics of the remaining moments couple to gradients of all the others. This 
shows that a judicious choice of moments (including negative ones where appropriate) is central to the simplification 
of the dynamics. 

The study of the linear dynamics shows also the appearance of a fast a a slow time scales (as was the case for 
chemical polydispersity), although in this case the role of the fast and slow modes are inverted. A detailed analysis of 
the modes deserves further work to elucidate the structure of the slow time scale. 
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APPENDIX A: EIGENVALUES FOR CHEMICAL POLYDISPERSITY 

For a chemically polydisperse solution in the presence of a passive solvent, the two relevant eigenvalues in the 
linearized regime are obtained from eq ^l[ The only coupled moments are mo and mi. Therefore, the relevant linear 
system of equations is then defined by the matrix, 

/ -q 2 D {1 - mo [1 - (1 - mo)(r - 1) - r [1 + m (r - 1)] W {q)]} 2 X D N p mi(l - mo) {1 + m (r - 1)} W (g)V \ 



\ q 2 D mi {1 - (1 - mo)(r - 1) - r [1 + m (r - 1)] W (q)} -q 2 D {l - 2 X N P [rh 2 - m\(l - (r - 1)(1 - mo))] W (q) 2 } J 

Using Mathematica, we get its corresponding eigenvalues 

D a 2 

A = -^-{-l-(l+m (r-l)) 2 + 2 X [m 2 +m2(-l + (l-m )( r -l))] +A} 
+ {4x(-TO2 + m\(2 + mo(— 1 + r) - r)) + m (l + mo(— 1 + r))r} 



(Al) 



48 
48A 

2x(-4 + 5r + m (-l + r)(-2 + 3r))(rh rh 2 + fh\{-2 + m (2 - m + (-1 + m )r)))} (A2) 



{ 77^(1 + m (-l + r))(2 + m (-l + r))(-l + r)r + 8x 2 (m 2 + mf (-2 + m + r - m r)) 2 



Ai = {— 1 — (1 + m Q (r - l)) 2 + 2x [m\(r - 1)(1 - m ) + m 2 - rnf] - A} 
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D q 4 



-+ {4x(-m 2 + m\{2 + m (-l + r) - r)) + m (l + m (-l + r))r} 

{(mg(l + m (-l + r))(2 + m (-l + r))(-l + r)r) - 8x 2 (m 2 + m?(-2 + m + r - m r)) 2 



48 
A)<7 4 



48A 

+ -2x(-4 + 5r + m (-l + r)(-2 + 3r))(m m 2 + m\{-2 + m (2 - m + (-1 + m )r)))} (A3) 
where the quantity A can be expressed as 

A = {(l + (l + m (r-l)) 2 + 2x[-m 2 +ml-m\{r-l)(l-m Q )]) 2 



+ 4(1 + m (r - 1))(-1 - m (r - 1) + 2 X [m 2 (l + m (r - 1)) - m 2 r] )} 



1/2 



(A4) 



APPENDIX B: GENERAL EXPRESSION FOR FLUXES IN LENGTH-POLYDISPERSE CASE 



In the main text we have given only the expressions for the diffusive fluxes in the linearized case because the generic 
expressions become quite lengthy, and display them here instead. The generic expressions for the diffusive flux of 
species i can be expressed as 



1 

Nf 

+ [m-i 
+ Pi 



1 / 1 - m 



k v2 + 2i y2wo 



Vm_ 



Wo 
1 



1 



- - -(mo - ^V 2 m x + i^V 2 ) - m_r + lv 2 m - ^V 2 



V72 , 1 ~ m 

2i v mo + ^r 



1 vt2 1 vr2 

m ° ~ 24 mi 24 TOl 



1 



247V, 



-V 2 



m_i + 



l-m \ 1 2 



24 



1 / 1 m 1 _ 2 

^U- + ^- m - 1+ 24 V m ° 



Vm 



while the corresponding fluxes for the moments can be expressed as 

m„(l - mo) 



Vm„_ 2 + m„ Vm_ 2 + 
2X 



m„_i - m n m_i 
m n m 



+ ^™«(V 2 m - V 2 ) 



Vm_ 



+ m„m_i + 



(1 - mo) ^m„_i - m„m_i + 
m„(l - m ) 



[V 2 m„ + m„((l - m )V 2 mo - m V 2 ) 



1 

+ 24 



N s 

m„_i + m„TO_i 



m„_ij (miV - V mi) 
m„(l - m )^ ^2 



m„_i - m„(m_i - m /N s ) 



JV, 



247V, 



V m > Vm 



(Bl) 



(B2) 



APPENDIX C: DIFFUSIVE FLUXES AFTER DYNAMICAL PROJECTION 



We provide here the expression for the diffusive fluxes of the moments m_2, m_i and mo, when we have derived a 
closed set of equation for their dynamics. These expressions were used to obtain numerically the three corresponding 
eigenmodes, and can be written down as 



1 



j-2 = (-W 2 742 + m (iV)(iV- 1 ))^Vm_ 2 

+ [-731 (N) 3 + m {N- 2 )(N) 2 - »no<JV- 1 )<JV)(r(l - m ) + m ) 



Vm_! 



+ [- (N) 4 T40 + rh (2 X (N)(l - m ) - r)((7V- 2 )(iV> 2 + m (r - l)^" 1 )^))] 



Vm 



j-i = {-(N) T 32 + ™ )^yVm_ 2 + (-W 2 73i + mo(JV _1 )W - m ((l - m )r + m ))^Vm_ 1 
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+ [-(N) 3 T30 + 2x(N)m (l - ^((iVr 1 ) W + m (r - 1)) - ^(^-^(JV) + rh (r - 1))] 



1 



(7V)3 



Vtoq 



(1 - mo)Vm_2 - m (l 



m )(r 



1) 



Vm_i 



Jo 



(TV) 



The coefficients 717 with 



+ m [2x(N)(l-m )(l + m (r 




the moments m„4 and to_3 into the subspace spanned by the moments to_2, and uiq. The values of these 

coefficients are functions of the subspace chosen. 
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Figure captions 

1. Figure 1: Wavevector dependence of the unstable eigenmode XiN p /D for an exponential parent p n ~ exp(a/„) 

— 1/2 

at x = 1.2\s- The wave vector q is measured in units of N p . In all the plots, we show the eigenmodcs for 
four different size ratios r = 10, 50, 100, 300 and 1200. In the different plots we change the degree of asymmetry 
of the chemical distribution (increasing with increasing a), and compare dense and dilute suspensions, a) a = 2, 
too = 0.15; b) a = 5, m = 0.15; c) a = 10, to = 0.15 ; d) a = 10, rh — 0.8. 

2. Figure 2: Wavevector dependence of the stable eigenmode XqN p /Dq for an exponential parent p n ~ exp(a/ n ) 

at high dilution m = 0.15 and \ — l-2Xs- The wave vector g is measured in units of N p 1 . a) a = 2; b) 
a =10. 

3. Figure 3: Quenched (long-dashed) and annealed (continuous) spinodals for a chemically 

polydisperse mixture with an exponential parent shape p n ~ exp(a/ n ) as a function of monomer concentration 
for various degrees of asymmetry a and size ratio r. The interaction strength \ is normalized by its value at 
the critical point, Xc for the corresponding set of parameters r and a. a) a = 1, r = 10; b) a = 5, r = 10; c) 
a = 20, r = 10;d) a = 20, r = 1000. 
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Figure 2b 
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Figure 3 a 
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Figure 3b 
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Figure 3 c 
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Figure 3d 




